Anharmonic oscillation effect on the Davydov-Scott monomer in thermal bath 
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The dynamics of Davydov-Scott monomer in a thermal bath with higher order amide-site's dis- 
placement leads to anharmonic oscillation effect is investigated using full-quantum approach and 
the Lindblad formulation of master equation. The specific heat is calculated based on the thermo- 
dynamic partition function using the path integral method. The temperature dependence of the 
specific heat is studied. In the model the specific heat anomaly as pointed out in recent works by 
Ingold et.al. is also observed. However it is found that the anomaly occurs at high temperature 
region, and the anharmonic oscillation restores the positivity of specific heat. 

PACS numbers: 87.10.-e, 03.65. Yz 



c 

o 

(N 
> 

o 
o 

o 
o 



X 



I. INTRODUCTION 

The study of molecular biophysics began with Frolich's 
hypotheses which assumes that an excitation from an 
atomic vibration related to biological activity [1]. Based 
on this idea, Davydov has developed a quantum theory 
of protein to understand the mechanisms of energy trans- 
port in molecular protein, in particular alpha helix pro- 
tein. The mechanism is an excitation energy of an amidc- 
I is stabilized by its vibration in a combined excitation 
which propagate as a soliton [2, 3]. Most of studies in 
this field have been done at zero temperature, and little 
attention has been given at physiological temperatures. 

The studies of Davydov soliton in physiological tem- 
perature is realized by the Davydov soliton interact- 
ing with thermal bath. While zero temperature cal- 
culation allows the existence of a soliton states in the 
proteins, there is an important question whether these 
states are stable at biological temperatures [4, 5]. Be- 
cause the temperature effect may exchange heat energy 
with surrounding aqueous medium. The measurement of 
infrared absorption and Raman's scattering of an crys- 
tallineacetanilide (CHACON HC§H§) X at low tempera- 
ture showed a new band closing to the amide-I band [6, 7]. 
The result is interpreted as a signature of Davydov's soli- 
ton. The experiment using femtosecond IR spectroscopy 
realizing a band of the amide-I from acetanilide (ACN) 
and N-methylacetanide (NMA) shows the dependencies 
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of absorption spectrum on temperature. At high temper- 
ature, the absorption spectrum shifts to higher frequency 
[8] . Theoretical prediction using first order perturbation 
methods based on Davydov model gives soliton life time 
at the order of 0(1CU 12 ) s. This is very short time for 
biological processes in room temperature [9] . 

Using standard Davydov model, some numerical calcu- 
lations showed that soliton is stable at 310K [4]. The cal- 
culation based on trial function by Kapor et.al. showed 
that soliton is stable at 300K [10]. The result has also 
been confirmed in [11]. 

From those results, it is important to study Davydov 
model involving the contact with thermal bath. The be- 
havior of the system has attracted many interests in the 
last three decades. The interaction of a system with its 
environment is given by the dissipation effect in quantum 
system [12-14]. However, the dissipation effect leads to 
a serious problem for quantization procedure due to the 
broken Heisenberg's relation. The most appropriate the- 
ory to resolve this problem is the Linblad formulation of 
master equation [13, 14]. 

The first application of Linblad formulation of master 
equation to the protein model has been done by Cuevas 
et.al. [15]. They used Davydov-Scott monomer model 
and showed that at 10K the quantum effect of amide-I vi- 
bration can not be neglected. Also at room temperature 
the semi classical approach might be a good approxima- 
tion compared to the corresponding full quantum system. 
However the study was focused on the dynamical aspect 
of the system. Since the real world is affected by thermal 
fluctuation, it is important to study such system from 
statistical mechanics point of view. 

In this paper the thermodynamic properties of 
Davydov-Scott monomer is investigated using Lindblad 
formulation and the partition function is calculated us- 
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ing path integral method [16]. The path integral method 
is a powerful tool to investigate the properties of non- 
linear dynamical systems with retarded interactions [17]. 
In Sec. II the Hamiltonian of the system under consid- 
eration is described as a coupled harmonics oscillation 
i.e. the amide-I and amide-site oscillators. In addition 
to the original Davydov-Scott monomer, the higher or- 
der of amide-site displacement inducing the anharmonic 
oscillation effect of amide-site is also taken into account. 
This is motivated by the fact that the excitation and 
relaxation of collective modes of a protein are generally 
achieved via anharmonic interactions with other normal 
modes through energy exchange [18]. In Sec. Ill the par- 
tition function is calculated using path integral approach 
to obtain the thermodynamic specific heat of the system 
which is then discussed in Sec. IV. 



II. HAMILTONIAN OF DAVYDOV-SCOTT 
MONOMER WITH ANHARMONIC 
OSCILLATION EFFECT IN LINDBLAD 
FORMULATION OF MASTER EQUATION 

Considering the Davydov-Scott monomer, the excita- 
tion of the amide-I is described by the coordinate (x) 
and momentum (p) operators. On the other hand, the 
displacement and momentum operators of amide-site are 
expressed by Q and P. The hamiltonian for Davydov- 
Scott monomer can then be written in the full quantum 
approach as follow, 

n 2 1 IP 2 f)V 

where to is the intrinsic frequency of amide-I oscilation, 
X measures the coupling between the amide-I excita- 
tion and the amide site vibration, m (M) is the amide-I 
(amide-site) mass and 8 is the anharmonic coefficient. 
Here the amide-site potential can generally be expanded 
around the origin, 
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The first term in Eq. (2) is a constant and can be scaled 
out. On the other hand, the equilibrium condition im- 
plies that the term of Q vanishes due to — W(Q)\q=o = 
0. The third term is the harmonic oscillator one, while 
the remaining higher order terms represent the anhar- 
monic contributions. 

The original Davydov-Scott model assumes that the 
anharmonicity of amide-site is not important compared 
to the amide-I [19]. The nonlinear effects appear only 
through exciton and amide-site coupling and higher or- 
der potential of excitons. However, recent experiment 
in biopolymer showed that the anharmonicity of poly- 
mer amide-site might be important [18, 20, 21]. The 



solid state experiments showed that at high temperature 
(room temperature), specific heat is not a constant any- 
more, but it is increasing. This fact can only be explained 
by taking into account the anharmonic oscillator term 
[22]. 

In Eq. (2), the terms Q 3 , Q 5 , • • • represent the asymme- 
try of atom mutual repulsion, while the terms Q 4 , Q 6 , ■ ■ ■ 
represent the softening of vibration at large amplitudes 
[22]. Moreover, concerning the fact that the Davydov- 
Scott polymer has an identic monomer, it is plausible 
to assume that the asymmetry of atom mutual repulsion 
can be ignored. This yields, 



V{Q) = \kQ 2 



I^ 4 



(3) 



where k/2 = 1/2 (d 2 V / dQ 2 )\ Q=0 and 

l/24(d 4 V/dQ 4 )\ Q=0 = -A/4. Again, the second 
term induces the anharmonic oscillation effect in the 
system. 

Unfortunately the term leads to a severe problem of 
describing the damping in open quantum systems, and 
it has been discussed for a long time. One of the known 
models dealing with this problem is the one dimensional 
damped harmonic oscillator (known as the Caldirola- 
Kanai Hamiltonian) [23]. In this model the momentum 
and coordinate operators are multiplied by e~ 7 * and e 7 *, 
where 7 is a damping factor. 

In statistical mechanics, the behavior of an open sys- 
tem within system-plus-bath can be modeled by the den- 
sity matrix formalism p. The equation of density matrix 
with hamiltonian H and environment operator R satisfy 
particular master equation [12]. It is usually restricted 
to weak system-bath interaction. The density operator 
gives the probability for the expected outcomes of mea- 
surements on the system. However, this formulation does 
not preserve density operator properties, that is hermitic- 
ity, unit trace, and positivity. The open quantum system 
theory which preserves density matrix properties can be 
realized using Lindblad formulation. In this theory, the 
interaction between Hamiltonian and thermal bath is re- 
alized by introducing some operators called Lindblad op- 
erators. The operators obey the master equation [24], 
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where Lj are the Lindblad operators. This choice is not 
unique and not necessarily Hermitian. Since L must be 
the first order in Q and P [14]. The operators H and 
L denote the internal dynamics and environmental ef- 
fects of the system. Throughout the paper, the Lindblad 
operators are put, 



Li = VW+ 




(5) 
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where 7 is a damping parameter related to the intensity 
of thermal bath, 0, — k/M, ks is the Boltzmann con- 
stant, T is temperature and v — ( e hn l k BT _ -j^-i j g ^ e 
Bose-Einstein distribution. 

Substituting Eqs. (5) and (6) into Lindblad master 
equation in Eq. (4), 

-^[x,Pp + pP]-^[Q,[Q,p]] 
-^ [x ' [Q ' p]] - lk [x ' [x ' p]] ~ ^h [P > [P ' p]] ' (7) 

where <Vs are the coefficient in the Lindblad operators, 



where 
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<5 12 are the frictional damping rate, while £3,4,5,6 are the 
quantum mechanical diffusion coefficients [25]. This is 
the underlying model in the paper. 



III. PATH INTEGRAL CALCULATION OF 
PARTITION FUNCTION 

In order to solve the master equation in Eq. (7), the 
resolution of a set of differential equations among the 
matrix elements of density operator with respect to a 
specific basis is required. This is usually provided by 
the eigen states of the system hamiltonian [26]. But in 
the statistical mechanics it is not necessary to solve the 
master equation. Instead one can calculate the partition 
function Z using, for instance, path integral method. 



A. Partition function 

The density matrix can be obtained by performing a 
transformation t — > r = —it = fih with (3 = l/(fceT) 
[16]. We assume that the quantum mechanical diffusion 
is dominant than the frictional damping rate such that 
it can be ignored. Then, the Lindblad master equation 
of the unnormalized p can be rewritten as, 



dp 

dp 



Hp + Ap 



(14) 



2h XP+ 2h Pp - 



(15) 



It should be remarked here that, using this equation one 
can confirm the Linblad operators in Eqs. (5) and (6) 
lead to the right equilibrium. The proof is given in App. 
A. 

The partition function corresponding to the master 
equation is given by [16], 



V[x{t)]e-* s(x(t)) 



(16) 



where S(x(t)) = Jq R (T + V)dt' is the Euclidean action 
corresponding to the equation. 

Our interest is on the anharmonic amide-site interac- 
tions. Choosing the amide-site potential as Eq. (3) and 
taking into account only the first order of potential in the 
interaction, the action becomes [16], 

f T /l 1 1 

S X Q — / dt I -mi 2 + —muj 2 x 2 — —Sx 4 + x x Q 
Jo \2 2 2 



+ 1 M Q 2 + \nQ 2 - ^Q 2 - \\Q 4 
2 2 2 



(17) 



This yields, 

StO — 



dt [ \mx 2 + \kx 2 - ^-Sx 4 + xxQ 



+ l -MQ 2 + l -RQ 2 \\Q 4 



(18) 



where k = k — iS 3 , k = k — iS 5 , M = M — iSe, 
X = X ~ i/2^4 and k = muj 2 . The amide-site is as- 
sumed to be more rigid than the amide-I. So the quan- 
tum fluctuation is dominated by the amide-I to enable 
us to use the Gaussian approximation. Making use of 
the Gaussian approximation, only the classical path of Q 
contributes to the interaction term [27]. Therefore the 
partition function becomes, 



where 



Z — Z x Zq , 

Z x = J Vxer s -' h , 
Zq= J VQc- s ^/ h , 



(19) 

(20) 
(21) 



and the actions are, 



S x = j dt Qmi 2 + X -kx 2 - l -bx 4 + xxQ^j , (22) 
S Q = £ dt QmQ 2 + \kQ 2 - \\Q^j . (23) 
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B. Partition function for the amide-site 



where, 



Further, one should solve the partition function of the 
amide-site (Zq). We use Gaussian approximation to 
solve the partition function. Under this approximation 
the general path can be expressed in the usual way as 
Q = Q + Q, where Q is the classical path and Q is the 
quantum path [16, 27]. Expanding the action in Taylor 
series Sq becomes [27], 



1 



1 



(24) 



Since the classical path satisfies the variational principles, 
SSq = 0, and taking only the second order, 



Zq = j VQexp 
This yields, 
Z Q = c - s Q/ h [ VQcxp 

+2 



1 f , 1 

-- / dt— 

hi 2 



'd 2 C cV 



=■ QQ 



dQQ 



v dQ , 



Q 



where, 
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Calculating the second variation, and substituting the 
result into Zq in Eq. (26) one has, 



e *«' n Z Qo 



-s%/h 
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ir 

hJo 



1 -J2 1 _2 1—4 
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x / VQcxp 



1^ 



1 

h 



dt 



x ( t;MQ + ^-RQ 2 - lxQ 2 Q 
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The equation of motion for Q in Euclidean coordinate is 
given by, 



A 2 = ~Qt-±imQl 



(32) 
(33) 



Now the problem is turned into solving the prefactor 
in the path integral, 



1 



.1; ( l -MQ 
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Using semi-classical approximation it can be rewritten as 

[28], 
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and its second order variation is given by, 
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The result is. 
5 2 S A 
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271" h Jo 
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dt MQ + (-k + 3AQ 2 ) Q 



and by substituting Q in Eq. (30), 



(37) 



S 2 S, 



Qo 



SQ 2 



(3VMk 
2tt 



[Aisech 4 (Ofi/3) tanh(fifi/3) 



x (cosh(2fifi^) - 3) + A 2 tanh 3 (f>ft/3) 
+A 3 sech 2 (nfi^) tanh 3 (TWi/3) 
x(4 + cosh(2fifi^)] , 



(38) 



where, 





Ai = 


\mqI, 

— — kVlQ a 
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A 2 = 


where Q = k/M. In this paper the solution is taken to 




have the form of, 


A 3 = 





Q = Qq sech(Qt) 



(30) 



which leads to Q = (2M)/X. This choice is moti- 
vated by the fact that the Davydov model is the self- 
trapping of its energy, i.e. it should be localized. Sub- 
stituting this solution into classical action and using the 
identity 1 — scch 2 a; = tanh 2 x, one obtains, 

S Q ] = A l tanh (Qh/3) + A 2 tanh 3 (Qh/3) , (31) 



(39) 
(40) 
(41) 

Finally, the complete partition function for the amide- 
site is obtained, 

Z Q = y^J\ C^«) 1/4 {Aisech 4 (ITfi^)tanh(nft/3) 
x (cosh(2Q/i/3) - 3) - A 2 tanh 3 (ftft/3) 
+A 3 sech 2 (n^/3) tanh 3 (Qh/3) (4 + cosh(2n^/3))} 1/2 



xcxp 



-jiM tanh(Q/i/3) + A 2 tanh 3 (Qfi/3)) 



.(42) 
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Partition function for the amide-I 



For the action of amidc-I, S x , the partition function is, 



= J Vxexp — ^ J dt ^ 
—ox 



1,1-, — 

-mx z + -kx z + \xQ 



(43) 



Dividing x into classical path x and quantum path x, i.e. 
x = x + x, and again using the Gaussian approximation, 
the classical path is, 



(44) 



This can be solved by determining the classical path 
which is the solution of following equation, 



rax — kx + 2Sw = xQo sech(Qi) . 



(45) 



However this is hard to be solved analytically. Instead 
one can use the perturbation method, i.e. calculating the 
solutions order by order x — x° + ex 1 + Note that 
the inhomogencous term is assumed being generated from 
the leading order. Substituting this expansion into Eq. 
(45) up to the leading orders one has, 



mx° — kx° 



2Sx 



-03 







(46) 



for the lowest order and, 



mx 



kx v + GSx^x 1 = xQo sech(Qt) , (47) 



for the first order. Similar to the solution in the amide- 
site coordinate, since S < k, the solution for the zeroth 
order is, 



Xq sech(u;t) 



(48) 



where ZJ = k/m, Xq = after transforming t — > 

it. In this solution 8 must not be zero. For the first order, 
the solution is x 1 = x\ + x p and x\ satisfies, 

§h - uP-xX + 6—^o scch 2 (uJi)^ = . (49) 

Performing a transformation r = tanh(uJi) one gets the 
associated Legendre equation, 



_d_ 

d7 



(1_ T 2)M 

dr 



+ + 



1 



0, (50) 



where n = uj and 1(1 + 1) = 6S/(muj)X . The solution 
is x\ = xj lQ P[ l (ta,nh.(LJt)) where P"(t) is the associated 
Legendre function. Particularly its solution satisfies, 

mXp — muJ 2 Xp+6SXQsech 2 (uJt)xp = xQo sech(f2t) . (51) 



The solution of this equation can be written using Green 
function G(r, r') as follow, 



x p — 



J dr' G(t, r')sech 



^tanh-^r') 



and the Green function is governed by, 



_d_ 

dr 



(i _ r2) dG(r,r') 



dr 



1(1 + 1) 



m 



1-t 2 



-8(t-t') . 



, (52) 

G(t,t>) 
(53) 



The Green function is given by [29], G(t,t') — 
(-1)™P ; "(t < )Q^(t > ) with Qf (t) is the associated Leg- 
endre functions of the second kind. Its complete solution 
is, 



x = X sech(wi) + e [x^P? 1 (ta,nh(uit)) 

+ J dr'G(T, rOsech^ltanh- 1 ^') 



.(54) 



Substituting this result into Eq. (44) one obtains the 
classical action. On the other hand, the classical action 
up to the first order is given by, 



Ai tanh(cJfi^) + A 2 tanh d (uJh/3) 
+A 3 / dt scch(uJt) sech(Qi) 



(55) 



+s(F 1 +F 2 + F 3 + F 4 + F 5 +F 6 + F 7 + F s ) , 



where, 



. 1 _- T7 i S =A 
A 1 = - m ^X + -X , 
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A 3 = -^xXoQo , 



(56) 

(57) 
(58) 



and Pi's are given in App. B. 

Performing the same procedure as done in the previous 
subsection, one should consider the prefactor of Z x , that 
is, 
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where the second order variation is given by, 
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which yields, 

5 2 S»a 



Sx 2 



2tt 



dt 



mx~ — {mui z + 35x 2 )x* 



(62) 



Substituting Eq. (54) into Eq. (62) and keeping only the 
first order, 



Sx 2 



/3V~r. 



2tt 



[T 1 sech 4 (wfy3) tanh(um/3)(cosh(2CJ/i/3) - 3) 
+r 2 tanh 3 (uJ/i/3) 

+r 3 (4 + cosh(2w^))sech 2 (aJ^) tanh 3 (aJ^) 
+e(G 1 + G 2 + G 3 + G 4 + G 5 
+G 6 + G7 + G 8 + G 9 + Gio)] , . 



where. 



ri-- y w X 



o > 



1TILJ 3 Xq 



r 2 = 



SLdXn 



and Gj's are given in the App. C. 

Hence, the partition function Z x is obtained, 



(63) 

(64) 
(65) 
(66) 



2-r — 



fiy/mu) 2 
4tt 2 H 



1/2 



[-r x scch 4 (wft/3) tanh(wft/3)(cosh(2aJ^) - 3) 
+r 2 tanh 3 (w/i/3) 

+r 3 (4 + cosh(2ZJhf3)) sech 2 (w/i/3) tanh 3 (ZJh/3) 
+e(Gi + G 2 + G 3 + G 4 + G 5 
+G 6 + G 7 + G 8 + G 9 + G 10 )] 1/2 

i |-i [Ai tanh(wfi/3) + A 2 tanh 3 (uJfi.,3) 
,/3fc 

\ / dt sech(wi) sech(Sli) 



xcxp 



+e(F 1 +F 2 + F 3 + F 4 + F 5 + F 6 + F r + F 8 )}} . (67) 

Finally, the complete partition function density of 
Davydov-Scott monomer in thermal bath, 2 — 2 X Zn , is 
obtained from Eqs. (42) and (67). 



IV. THERMODYNAMIC AL PROPERTIES OF 
DAVYDOV-SCOTT MONOMER 

From physical point of view, the Davydov-Scott 
monomer is a harmonic oscillator coupled to a quantum 
excitation. Using Euler-Lagrange equation, one can de- 
rive the appropriate equation of motion (EOM) from the 



action in Eq. (17). Although solving of the EOM is very 
interesting and attractive, but it has a little physical sig- 
nificant due to unobservable individual molecular motion 
of the Davydov-Scott monomer. Therefore in this paper 
let us consider thermodynamic observable as specific heat 
[16], 



C = k B f3 



2 d 2 ln(Z) 

d/3 2 



(68) 



Some previous works considering similar model as the 
Hamiltonian in Eq. (1) to study the Davydov-Scott 
monomer in thermal bath, for example the semiclassi- 
cal approach in [15]. It has been argued that using Lind- 
blad formulation the semiclassical limit is a good approx- 
imation to the corresponding full quantum treatment at 
biological temperatures in the highly underdamped and 
harmonic limits. In the semiclassical approximation, the 
coupling between Davydov-Scott monomer with thermal 
bath is described by Langevin equation of the amide- 
site displacement characterized by 7 and f2 [5]. If the 
stochastic force represents the thermal bath is zero, the 
equation is reduced into the damped harmonic oscilla- 
tor. There are three regions regarding the values of 7 
and 17, that is 7 < 2S7 for the underdamped condition, 
7 = 2fi for the critical damped condition and 7 > 2£1 
for the overdampcd condition. In the Lindblad formula- 
tion the damping coefficient 7 represents the relaxation 
time due to interaction with the environment. The higher 
values of 7 corresponds to the shorter relaxation time. 
The previous work by Cuevas et.al. has established that 
the semiclassical approximation is equivalent to the full 
quantum approach (for biological temperature) as long 
as 7 <C 2f2. Otherwise, the oscilation frequency of the 
observable would be different [15]. 

In this paper, the analysis is done for the above three 
criterions. The values of parameters used throughout 
numerical calculation are M — 6.3 x 10~ 26 kg, m = 
7.3 xlCT 26 kg, k = 10 Nm- 1 and u = 1660 cm- 1 [15, 30], 
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FIG. 1. The temperature dependence of normalized specific 
heat for various values of 7 with A = 1. 
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FIG. 2. The temperature dependence of normalized specific 
heat for various values of A for the underdamped condition 
with 7 = O.lfi. 



FIG. 4. The temperature dependence of normalized specific 
heat for various values of A for the overdamped condition with 
7 = 2.5Q. 



while A = 9. The behavior of specific heat in term of 
temperature is shown in Fig. 1 for various damping pa- 
rameter 7, Figs. 2—4 for various strength of anharmonic 
oscillation in three damping conditions and Fig. 5 for 
various 7 in the underdamped case. The results are sim- 
ilar with the previous ones obtained in the calculation of 
a system with anharmonic amide-site [31]. 

In the present case the damping coefficients are repre- 
sented by the coefficients Si <~ Sq. In particular, Si ap- 
pears in the kinetic terms, and can then be interpreted as 
the 'effective mass' of amide-site vibration. Further, S3 
appears in the harmonic potential as the 'effective elastic 
constant'. Hence it can be argued that the environment 
effects to the amide-site vibration occur through the ki- 
netic term and the harmonic potential. The coefficient S4 
represents the strength of the interaction between amide- 
Land the system. 



From the figures, the specific heat asymptotically ap- 
proaches to zero at low temperature and to infinity at 
high temperature. Actually <j) 4 potential oscillator in- 
teracting with electron and also the anharmonic oscilla- 
tor give similar profiles [32]. Large environment effect 
causes the Davydov-Scott monomer to increase the en- 
ergy and then the temperature to achieve the equilibrium 
state. Correspondingly the vibration frequency is also af- 
fected by the damping coefficients since u> = 

\J{k + iSz)/(M + iS~2). These results indicate that the 
interaction between Davydov-Scott monomer and ther- 
mal bath depend on the strength of the coupling of sys- 
tem and environment. Recent study of the open quan- 
tum system also shows that the canonical equilibrium 
state of an open quantum system depends explicitly on 
the system-bath coupling strength [33, 34]. 




j , , I -0.2 1 1 1 1 ~ 1 1 

l4 06 08 ! 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

W ' ' Wo 



FIG. 3. The temperature dependence of normalized specific 
heat for various values of A for the critical damped condition 
with 7 = 2Q. 



FIG. 5. The temperature dependence of normalized specific 
heat for the underdamped condition with various values of 
7 < 2Q and A = 1. 
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In particular, from Fig. 2 the anomaly of specific heat 
that becomes negative for certain parameter sets at high 
temperature region is observed. The same phenomena 
have been pointed out by Ingold et.al. [33, 35]. In the 
current case, the anomaly especially appears for the un- 
derdamped condition as shown in Figs. 2 and 5. It is 
also found that the negative specific heat is restored at 
large A and 7, i.e. for large anharmonic oscillation and 
intensity of thermal bath. 

It should also be remarked that one cannot take A = 
(no oscilation effect) since all results are obtained from 
Eq. (30) which is a special case with A 7^ condition. 



petitif LIPI in fiscal year 2010 under Contract no. 
11.04/SK/KPPI/II/2010. FPZ is supported by Risct KK 
2010 Institut Tcknologi Bandung. 



Appendix A: The equilibrium state with the 
Lindblad operators in Eqs. (5) and (6) 

Let us investigate the equilibrium state in the present 
case. Substituting the amide-site potential in Eq. (3) and 
the leading interaction term in Eq. (3) into the Lindblad 
master equation in Eq. (14) yields, 



V. SUMMARY 

The interaction of Davydov-Scott monomer with ther- 
mal bath is investigated using the Lindblad formulation 
of master equation. In contrast with previous work by 
Cuevas et.al. [15], the anharmonic oscillation term of 
amide-site is taken into account. Adopting similar Lind- 
blad operators used in [15], the master equation of the 
system is obtained. Instead of solving the equation of mo- 
tion, the thermodynamic partition function and in par- 
ticular specific heat are calculated using the path integral 
methods. 

It is shown that the coupling with the environment 
contributes to the kinetic term, the harmonic potential of 
amide-site vibration and the anharmonic term of amidc- 
I. 

The anomaly of specific heat that becomes negative for 
certain parameter sets at high temperature region is ob- 
served as pointed out by Ingold et. al. in the case of pure 
open quantum systems [33, 35]. However, it is found 
that the negative specific heat is restored for large an- 
harmonic oscillation effect. In contrast to these results, 
Ingold et.al. have found that the anomaly occurs at low 
temperature region. This discrepancy can be explained 
as the consequences of different approaches adopted to 
model the interaction between the system and the ther- 
mal bath. Ingold et. al. has represented the interaction in 
a set of harmonic oscillators which becomes the coupled 
harmonic oscillator at classical approximation. In con- 
trary, in the present approach the thermal bath is repre- 
sented in a set of Lindblad operators which becomes the 
underdamped harmonic oscillator at classical approxima- 
tion. 

From the figures, it can in general be concluded that 
the anharmonic oscillation contributes constructively to 
the specific heat. 
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This can be rearranged to be, 
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On the other hand, the average value of any operator 
A at a thermal equilibrium is given by [16], 



or in the integral form, 

((A)) = 



Tr(e-( m A) 
Tr(e-f> H ) 



J p(x)A(x) dx 
J p(x) dx 



(A3) 



(A4) 



Before going further, it is more convenient to rewrite Eq. 
(A2) as follow, 



dp 

d/3 



j_(Pp 

2m dx 2 

1 d 2 p 
+ 2MdQ 2 



1 9 2 

—moo x p 



KQ 2 P+i\Q 4 P + x*Qp (A5) 



-Sx p 
1 



m + iS 5 /uj 2 , M = Mh/(h + 1S4M), k = k + 
- x + iSi/{2h). It is clear that the Lindblad 



where m ■ 
183 and x 

operators shift all of the parameters, while contribute 
nothing to the fourth power terms, i.e. the third and 
sixth terms. Therefore, for the sake of simplicity one 
can ignore the fourth power terms when investigating the 
thermal equilibrium condition. 
This fact leads to, 

dp 1 d 2 p 1_ 2 2 
-dp = ^dx^ + 2 mLJXP 



1 d 2 p 1 ^ 2 

+ mW 2 + r Q p + * xQp - 



(A6) 
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Making use of the Gaussian approximation as before, 
Q = Q + Q [16], and assuming that only the classical 
part of amide-I contributes to the interaction, one imme- 
diately obtains, 



dp 1 d 2 p : l_ 2 ., 



d/3 2m dx 2 2 

1 8 2 p 1 



mu x p + xQ x P 



+ 



2M dQ 2 2 



(A7) 



This equation can be splitted into two equations belong- 
ing to the regular, 



dp _ 1 d 2 p 1 " 2 
and the driven oscillator harmonics, 



(A8) 



dp 1 d 2 p 1_ 2 2 - /4 . 

-dj3^2^dx^ + 2 mUJX P + XQXP - (A9) 

The solution for Eq.(A8) is [16], 



MQ, 



2vrftsinh(nO/3) 



exp 



~ tanh ( Imp ] (f 



(A10) 

with £1 = \Jk/M. Subsequently, the thermal equilibrium 

for Q 2 can easily be calculated using Gaussian integral 
to get, 



((Q 2 )) - JL= coth ( hn/3 
U " 2Mfl \2 



The internal energy is given by, 
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e = — coth -mib = 

2 \2 ' J 2 i- e -hcip 

_ Ml Mle m P 
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(All) 
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Meanwhile, the oscillator harmonic with nq amide-site 
has the energy E = Ml{\ + uq). Hence, the number of 
quanta for amide-site at thermal equilibrium becomes, 



-ftfi/3 



(A13) 



as expected. Particularly, the case of 63 = S4 = repro- 
duces the oscillator harmonic at equilibrium without any 
environmental effects. 

Following the same procedure, one can obtain the ther- 
mal equilibrium condition for amide-I. Under the initial 
condition p(0) = 5(x — x'), the solution for Eq. (A9) is 
[16], 



p(x) = 



moj 



2-kH sinh(fto;/3) 



^(-Ax 2 +Bx+C) 



(A14) 



where, 



A = -^ muj tanh (^fuv(3 



(A15) 
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Then the thermal equilibrium for x 2 is, 

jc(- Ax2+Bx+c ')x 2 dx 



(A17) 
(A18) 

(A19) 



(A20) 



This integral is well known, and can be calculated by 
performing the transformation, x — B/2A and £ = x — x, 
and defining g(x) — B 2 /AA + C as well. These yield, 



((x 2 ))=c 



I 



e~ A i 2 ed(i + 2x- J e- A t 2 £dZ 



+x 2 I e"^ 



dt 
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(A21) 



Using the Gaussian integral, i.e. J e A ^ 2 £ 2 d£ = 

0F/(2A 3 / 2 ), Jc- A « 2 £d£ - and fe^d^ = y/H/y/A, 
the solution is, 



_ V¥/(2A 3 ' /2 ) + x 2 v / ^/v / I 



J_ ff 2 
2A + il 2 ■ 



(A22) 



Substituting Eqs. (A15) and (A16) yields, 



{{x )) = -z coth 

u " 2mw 



((xQ)) , (A23) 



where, 



((xQ)) 



. coth 2 -huj/3 
sinh(/zw/3) J V 2 

x [(A 1+ A 2 )(c^-1)] 2 , (A24) 



represents the coupling effect between amide-I and 
amide-site. The internal energy is given by, 



E=^- coth ) - 777^'" ( ,<): 
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Again, concerning that E = hw(l/2 + n x ), the number 
of quanta for amide-I at equilibrium becomes, 



n x = 



e haj/3 _l ft 



muJ 1 1 ~7\\ \ 

+ —((xQ)) 



(A26) 



The case of 5 3 = <S 4 = reproduces the number of quanta 
for amide-I at thermal equilibrium without any environ- 
mental effects. 

These results confirm that the Lindblad operators de- 
fined in Eqs. (5) and (6) lead to the right equilibrium as 
expected. 



Appendix B: The coefficients in Eq. (56) 
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Appendix C: The coefficients in Eq. (63) 
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